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Abstract. We present a written version of four lectures given at the NATO Advanced Study 
Institute on "QCD Perspectives on Hot and Dense Matter" in Cargese, Corsica during August, 
2001. Over the last year the first exciting results from the Relativistic Heavy Ion Collider 
(RHIC) and the four experiments BRAHMS, PHENIX, PHOBOS, and STAR have been pre- 
sented. In these lectures we review the state of RHIC and the experiments and the most exciting 
current results from Run I which took place in 2000. A complete review is not possible yet with 
many key results still preliminary or to be measured in Run II, which is currently underway, 
and thus the emphasis will be on the approach experimentalists have taken to address the 
fundamental physics issues of the field. We have not attempted to update the RHIC results for 
this proceedings, but rather present it as a snapshot of what was discussed in the workshop. The 
field is developing very quickly, and benefits greatly from contact and discussions between the 
different approaches of experimentalists and theorists. 
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1. Introduction 



A good way to proceed in understanding the construction, operation, and 
physics output from current high-energy heavy ion experiments is to start 
with the original physics motivations for the experimental program. The Rel- 
ativistic Heavy Ion Collider project was started in the 1980's with a list of 
experimental observables for characterizing hot and dense quark and glu- 
onic matter and the expected restoration of approximate chiral symmetry and 
screening of the long range confining potential of QCD. We give a brief and 
certainly not complete review of these physics signals. Then we discuss some 
of the cost, schedule and technology constraints that impacted the design and 
construction of the RHIC experiments. Finally we present a sample of the first 
results from the RHIC experimental program from data taken during Run I in 
the summer of 2000. 



2. Physics Motivations 

There are three main categories of observables that were originally proposed 
to study the matter produced in relativistic heavy ion collisions: (1) Thermo- 
dynamic properties of the system and indications of a possible first or second 
order phase transition between hadronic matter and quark-gluon plasma, (2) 
Signatures for the restoration of approximate chiral symmetry transition, and 
(3) Signatures for deconfinement and the screening of the long range con- 
fining potential between color charges. In the following list, there will be 
specific channels that need to be observed that essentially specify the types 
of detectors and experiments necessary. 

2.1. Thermodynamic Properties 

One of the most important question in the physics of heavy-ion collisions is 
thermalization. We want to describe the system in terms of a few thermody- 
namic properties, otherwise it is not possible to discuss an equation-of-state 
and a true order to any associated phase transitions. The use of thermody- 
namic concepts to multi-particle production has a long history. One of the 
first to apply them in elementary collisions was Hagedorn in the early 1960's 
[1]. The concept of a temperature applies strictly speaking only to systems in 
at least local thermal equilibrium. Thermalization is normally only thought to 
occur in the transverse degrees of freedom as reflected in the Lorentz invariant 
distributions of particles. The measured hadron spectra contain two pieces 
of information: (i) their normalization, i.e. their yields ratios, provide the 
chemical composition of the fireball at the chemical freeze-out point where 
the hadron abundances freeze-out and ( ii) their transverse momentum spectra 
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which provides information about thermalization of the momentum distribu- 
tions and collective flow. The latter is caused by thermodynamic pressure and 
reflects the integrated equation of state of the fireball matter. It is obvious that 
the observed particle spectra do not reflect earlier conditions, i.e. the hot and 
dense deconfined phase, where chemical and thermal equilibrium may have 
been established, since re-scattering erases most traces from the dense phase. 
Nevertheless, those which are accumulative during the expansion, such as 
flow, remain. Only direct photons, either real photons or virtual photons that 
split into lepton pairs, escape the system without re-scattering. Thus these 
electromagnetic probes yield information on the earliest thermodynamic state 
which may be dominated by intense quark-quark scattering. 

The assumption of a locally fhermalized source in chemical equilib- 
rium can be tested by using statistical thermal models to describe the ra- 
tios of various emitted particles. This yields a baryon chemical potential 
fiB, a strangeness saturation factor 7 S , and the temperature T c h at chemical 
freeze-out. 

So far these models are very successful in describing particle ratios at 
SPS [2] and now also at RHIC [3]. At RHIC the derived chemical freeze- 
out temperature are found to be around 175 MeV (165 MeV at SPS) and a 
baryon chemical potential of around 45 MeV (270 MeV at SPS). It should 
be stressed that these models assume thermal equilibration but their success 
together with the large collective flow (radial and elliptic) measured at RHIC 
is a strong hint that this picture indeed applies. 

2.2. Chiral Symmetry Restoration 

The phase transition to a quark-gluon plasma is expected to be associated with 
a strong change in the chiral condensate, often referred to as the restoration of 
approximate chiral symmetry in relativistic heavy ion collisions is discussed 
here. Note that although many in the field refer to the restoration of chiral 
symmetry, the system always breaks chiral symmetry at a small scale due to 
the non-zero neutral current masses of the up and down quarks. 

There are multiple signatures of this transition, including disoriented 
chiral condensates (DCC), strangeness enhancement, and many others. How- 
ever, the most promising signature is the in medium modification to the mass 
and width of the low mass vector mesons. Nature has provided an excellent 
set of probes in the various low mass vector meson states (p, lo, 0) whose mass 
poles and spectra are dynamically determined via the collisions of hadrons 
or partons. If the hot and dense state produced in heavy ion collisions is 
composed of nearly massless partons, the /o(770) meson mass distribution is 
expected to broaden significantly and shift to lower values of invariant mass. 
The p meson has a lifetime that is r « 1 fm/c and the plasma state created in 
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RHIC collisions has a lifetime of order 10 fm/c. Thus, the p meson is created 
and decays many times during the entire time evolution of the collision. 

The p has a dominant decay (nearly 100%) into two pions. However, in 
this decay mode if the pions suffer re-scattering with other hadrons after the 
decay, one cannot experimentally reconstruct the p meson and information is 
lost. Given the dense, either partonic or hadronic environment, the probability 
for pion re-scattering is very large, unless the p is created and decays at the 
latest stages of the time evolution. This time period is often referred to as 
thermal freeze-out, when elastic collisions cease. Thus a measurement of the 
p as reconstructed via its pion decay channel gives interesting information on 
the final hadronic stage, but not on the dense phase where chiral symmetry 
may be restored. 

There is an additional decay channel into electron pairs and muon pairs, 
though with small branching ratios of 4.5 • 10~ 5 and 4.6 • 10 -5 , respectively. 
Since the leptons do not interact strongly, after they are produced, they exit 
the dense system essentially unaffected carrying out crucial information from 
the core of the system. There is a good analogy in understanding the processes 
in the center of the sun via neutrino emission, since only neutrinos have a 
small enough interaction cross-section to pass out of the sun's core largely 
unaffected. In the case of neutrinos, the more interesting physics of possible 
neutrino oscillations complicates matters, but that is not a concern in the case 
of our electron and muons measurements. 

One other additional point of interest is that the apparent branching ratio 
of the p into pions and leptons should be modified as observed by experiment. 
The p mesons that decay in medium into pions are not reconstructed, but the 
ones decaying into electrons are. If one can reconstruct the p in both channels 
one can gauge the number of lifetimes of the p the dense medium survives 
for. 

The lifetime of the p meson in the rest frame of the plasma depends 
on its gamma boost in this frame, and thus to study the earliest stages, a 
measurement of low transverse momentum p mesons is desirable. If the p 
decays at rest in the plasma frame, the maximum transverse momentum for 
the electron or positron is p± ps 385 MeV/c. These electrons are considered 
low p± and present an experimental challenge to measure for two reasons. 
First, there are a large number of low momentum charged pions created in 
these collisions, that results in a charged pion to electron ratio in this p± range 
of 1000 : 1. Thus one needs detectors that can cleanly identify electrons with 
good momentum resolution, while rejecting the copiously produced pions. 
The second challenge is that most of the electrons come from pion Dalitz 
decays (ir° — ► e + e~7), r] Dalitz decays (77 — > e + e _ 7), and conversions of 
photons (7 — > e + e~) mostly resulting from 7r° decays. Ideally one wants 
to reject these other electrons to enhance the signal contribution from the 
low mass vector mesons. Conversions are reduced by reducing the amount of 
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material in the path of produced photons. This restriction is often at odds with 
the desire to have substantial inner tracking detectors and these needs must 
be balanced. 

The 0(1020) meson spectral function is also sensitive to in-medium 
chiral symmetry restoration; however its substantially longer lifetime r « 
40 fm/c means that most <p decays occur outside the medium. However, 
measuring the <p in its various decay modes (kaon pairs, electron pairs, muon 
pairs) remains an interesting signal at low transverse momentum. 

These low mass vector mesons also decay into muons pairs. These muons 
have low momentum and are an real experimental challenge to measure as 
detailed later in these proceedings. 

2.3. Deconfinement 

There are many signatures that result from the deconfinement of color charges 
over an extended volume, often referred to as the quark-gluon plasma. Two 
are detailed in this proceedings: (1) Suppression of heavy quarkonium states 
and (2) Partem energy loss via gluon emission, also referred to as jet quench- 
ing. 

2.3.1. Quarkonium Suppression 

The suppression of heavy quarkonium states was originally proposed by Mat- 
sui and Satz [4] in the late 1980's as a signature for color deconfinement. 
The Debye screening in a QED plasma is a reasonable analogue for the 
scenario in our QCD plasma. A charm-anticharm (cc) quark pair produced 
via gluon fusion in the initial phase of the heavy ion collision can form a J/ip 
if the pair has low relative momentum. The total production of such states in 
proton-proton collisions relative to the total charm production is less than a 
few percent. If the Debye screening length is of order the same size as the 
quarkonium state, then the pair is screened. The charm and anticharm quark 
scatter away from each other and, eventually at the hadronization point, pair 
with surrounding light quarks and antiquarks to form D mesons. This color 
screening is displayed in recent lattice QCD calculations described at this 
workshop in terms of a modification in the linear rise at large distances of the 
QCD potential. This change in the QCD potential as a function of temperature 
is shown in Fig. 1. 

There are a variety of heavy vector mesons with a large range in binding 
energy (and associated hadronic size). The J/tp, Xc, and tp' have binding 
energies of 0.64, 0.20, and 0.05 GeV respectively. The T(ls), Xfc> T(2s), 
X' b , and T(3s) have binding energies of 1.10, 0.67, 0.54, 0.31, 0.20 GeV re- 
spectively. Since the suppression of these states is determined by the relative 
plasma temperature and the binding energy (or by the quarkonium hadronic 
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Figure 1. The QCD potential between heavy quarks in shown as a function of ry/a 



size and the Debye screening length), measuring the sequential disappearance 
of these states acts as a QCD thermometer. 

The J/-0(ls) state decays into almost anything hadronic with a large 
branching ratio of 87.7%. However, the experimentally accessible decay chan- 
nels are 5.93% to e + e~ and 5.88% to Similar decay channels are 
available for the tp' and the T states. The accessible decay channel of the \c 
for heavy ion experiments is Xc —> 7 + J with a branching ratio 6.6 • 10 -3 . 
Since the decay 7 is quite soft, the \c represents an experimental challenge. 

Similar to the low mass vector mesons there is interest in these states at 
low transverse momentum, where they reside in the plasma state longer. At 
rest the J/ijj decays into electrons or muons with a characteristic p± as 1.5 
GeV. A rough rate estimate (good to a factor of 2-3) is that the production of 
J/ip is approximately 1 • 10~ 4 per proton-proton collision. In a central Au- 
Au collision, there are of order 800 binary collisions, yielding a J/ip rate of 
8 • 10 -2 . The branching ratio to electrons is 5.9% and a typical experimental 
acceptance is 1%, yielding 4 • 10~ 5 J/ip per Au-Au central collision, and that 
is assuming no anomalous suppression! Hence, one requires a detector that 
measures either electron or muon pairs with a high efficiency and a trigger and 
data acquisition system capable of sampling events at the full RHIC design 
luminosity. In particular, if one wants to bin the data in terms of xp, p±, and 
collision centrality, large statistics are a requirement. 

There is a recent proposal for J ftp enhancement. This scenario assumes 
copious charm prodcution, and then at the hadronization stage some charm 
and anticharm quarks may be close in both momentum and configuration 
space. These cc pairs may coalesce to form J/ip. This late stage production 
would potentially mask any suppression in the early stages from color screen- 
ing. There is an easy test of this theory. When RHIC runs at lower energies, 
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for example yfs = 60 GeV, instead of the maximum energy of yfs = 200 
GeV, the charm production is lower by a factor of approximately three and the 
effect of recombination should be reduced substantially. In addition, we ex- 
pect different p± dependence of J/ ip production from original hard processes 
compared with late stage cc coalescence. 

2.3.2. Parton Energy Loss 

An ideal experiment would be to contain the quark-gluon plasma, and send 
well calibrated probes through it, and measure the resulting transparency or 
opacity of the system. There is no experimental way of aiming a third beam 
of particle at the collision. Therefore any probes of the system must be gen- 
erated in the collision itself. These probes must have calculable production 
rates in order to be considered calibrated. An excellent example of such a 
probe is a hard scattered parton. A parton traversing a color confined medium 
of hadrons sees a relatively transparent system. However, a parton passing 
through a hot colored deconfined medium will lose substantial energy via 
gluon radiation[5, 6]. 

The source of these partons is from hard scattering processes producing 
back-to-back parton jets. In a deconfined medium the parton will lose energy 
before escaping the system and fragmenting into a jet cone of hadrons. The 
total energy of the initial parton jet is conserved since eventually the radiated 
gluons will also hadronize. It is likely that the radiated gluons will have 
a larger angular dispersion than the normally measured jet cone. Thus one 
might be able to measure a modification in the apparent jet shape. When the 
parton fragments into hadrons it has less energy, and hence the fragmentation 
will result in a much reduced energy for the leading hadron. A measurement 
of high transverse momentum hadrons (tt°, 7r + /~, K + /~ , is a strong 

indicator of the opacity of the medium. An exciting additional observable was 
mentioned at the workshop in the context of the high p± spectra of charm D 
mesons. There is a reduction in the induced gluon radiation for charm and 
bottom quarks relative to light quarks due to their slower velocity through the 
medium. 



3. The RHIC Complex 

The scope of the RHIC program is to operate a colliding beam facility which 
allows studies of phenomena in ultra-relativistic heavy-ion collisions and in 
collisions of polarized protons. The collider is located in the northwest section 
of the Brookhaven National Laboratory (BNL) in Upton, New York. Its con- 
struction began in 1991 and the completion of the complex was accomplished 
in Spring 2000. 
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The collider, which consists of two concentric rings of 1740 super-con- 
ducting magnets, was constructed in an already existing ring tunnel of ~ 3.8 
km circumference. This tunnel was originally constructed for the proposed 
ISABELLE project. It offers an extraordinary combination of energy, lumi- 
nosity and polarization. A schematic diagram of the whole RHIC complex, 
including the various facilities used to produce and pre-accelerate the beams 
of particles is displayed in Figure 2. 

RHIC is able to accelerate and store counter-rotating beams of ions rang- 
ing from those of gold to protons at the top energy of 100 GeV/nucleon for 
gold and 250 GeV for protons. The stored beam lifetime for gold in the energy 
range of 30 to 100 GeV/nucleon is expected to be approximately 10 hours. 
The major performance parameters are summarized in Figure 3. 




Figure 2. The Relativistic Heavy Ion Collider (RHIC) accelerator complex at Brookhaven 
National Laboratory. Nuclear beams are accelerated from the tandem Van de Graaff, through 
the transfer line into the AGS Booster and AGS prior to injection into RHIC. Details of the 
characteristics of proton and Au beams are also indicated after acceleration in each phase. 

The layout of the tunnel and the magnet configuration allow the two 
rings to intersect at six locations along their circumference. The top kinetic 
energy is 100+100 GeV/nucleon for gold ions. The operational momentum 
increases with the charge-to-mass ratio, resulting in kinetic energy of 125 
GeV/nucleon for lighter ions and 250 GeV for protons. The collider is able 
to operate a wide range from injection to top energies. The collider is de- 
signed for a Au-Au luminosity of about 2 • 10 26 cm~ 2 s _1 at top energy. This 
design corresponds to approximately 1400 Au-Au minimum bias collisions 
per second. The luminosity is energy dependent and decreases approximately 
proportionally as the operating energy decreases. For lighter ions it is signif- 
icantly higher reaching ~ 1 • 10 31 cm -2 s" 1 for pp collisions. The collider 
allows collisions of beams of equal ion species all the way down to pp and of 
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Figure 3. RHIC performance parameters. 

unequal species such as protons on gold ions. Another unique aspect of RHIC 
is the ability to collide beams of polarized protons (70-80%) which allows the 
measurement of the spin structure functions for the sea quarks and gluons. 

The first physics run at the Relativistic Heavy Ion Collider (RHIC) took 
place in the Summer of 2000. For this run beam energies were kept to a mod- 
erate 65 A GeV. RHIC attained its goal of ten percent of design luminosity by 
the end of its first run at the collision center-of-mass energy of ^/s NN = 130 
GeV. 

In the following we describe briefly the various facilities, depicted in 
Figure 2, that are part of the large RHIC complex: 

Tandem Van de Graaff Completed in 1970, the Tandem Van de Graaff fa- 
cility was for many years the world's largest electrostatic accelerator 
facility. It can provide beams of more than 40 different types of ions 
ranging from hydrogen to uranium. The facility consists of two 15 MV 
electrostatic accelerators, each about 24 meters long, aligned end-to- 
end.In the Tandem the atoms are stripped of some of their electrons (e.g. 
Au to Q = +32) and accelerated to a kinetic energy of 1 MeV/nucleon. 

Heavy Ion Transfer Line (HITL) To study heavy ion collisions at high en- 
ergies, a 700 meter-long tunnel and beam transport system called the 
Heavy Ion Transport Line were completed in 1986, allowing the delivery 
of heavy ions from the Tandem to the Booster for further acceleration. 
The HITL makes it possible for the Tandem to serve as the Relativistic 
Heavy Ion Collider's ions source. 
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Linear Accelerator (Linac) For the study of pp or pA collisions at the ex- 
periments, energetic protons are supplied by an Linear Accelerator (Linac) 
The Brookhaven Linear Accelerator was designed and built in the late 
1960's as a major upgrade to the Alternating Gradient Synchroton (AGS) 
complex. The basic components of the Linac include ion sources, a 
radiofrequency quadrapole pre-injector, and nine accelerator radiofre- 
quency cavities spanning the length of a 150 m tunnel. The Linac is 
capable of producing up to a 35 milliampere proton beam at energies up 
to 200 MeV for injection into the AGS Booster. 

Booster The Alternating Gradient Synchrotron Booster is less than one quar- 
ter the size of the AGS. It is used to preaccelerate particles entering 
the AGS ring and plays an important role in the operation of the Re- 
latavistic Heavy Ion Collider (RHIC) by accepting heavy ions from the 
Tandem Van de Graaff facility via the Heavy Ion Transfer Line (HITL) 
and protons from the Linac. It then feeds them to the AGS for further 
acceleration and delivery to RHIC. After the installation of the HITL in 
1986, the AGS was capable of accelerating ions up to silicon with its 
atomic mass of 28. However, due to its superior vacuum, the Booster 
makes it possible for the AGS to accelerate and deliver heavy ions up to 
gold with its atomic mass of 197. 

AGS Since 1960, the Alternating Gradient Synchrotron (AGS) has been one 
of the world's premiere particle accelerators and played a major role 
in the study of relativistic heavy ion collisions in the last decade. The 
AGS name is derived from the concept of alternating gradient focus- 
ing, in which the field gradients of the accelerator's 240 magnets are 
successively alternated inward and outward, permitting particles to be 
propelled and focused in both the horizontal and vertical plane at the 
same time. Among its other duties, the AGS is now used as an in- 
jector for the Relativistic Heavy Ion Collider. For RHIC operation the 
fully stripped ions are accelerated in the AGS to 9 GeV/nucleon before 
ejection. 

ATR The AGS sends the ions (or protons) down another beamline called 
the AGS-to-RHIC Transfer Line (ATR). At the end of this line, there's 
a "fork in the road", where sorting magnets separate the ion bunches. 
From here, the counter-rotating beams circulate in the RHIC where they 
are collided at one of four intersecting points. 
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4. Experimental Program 

4.1. Letters of Intent 

In July 1991 there were a set of experimental Letters of Intent that were put 
forward to an advisory committee. The proposals are listed below, including 
the lead institution and in parenthesis the physics observable focus. 

1. LBL-TPC (inclusive charged hadrons) 

2. BNL-TPC (inclusive charged hadrons) 

3. TOYKO-TALES (electron pairs, hadrons) 

4. SUNY-SB (direct photons) 

5. Columbia-OASIS (electron pairs, hadrons, high p±) 

6. ORNL Di-Muon (muon pairs) 

7. BNL Forward Angle Spectrometer (hadrons at large rapidity) 

8. MIT MARS (hadrons and particle correlations) 

The experiments span the range of hadronic, leptonic and photonic capa- 
bilities to cover the broad spectrum of physics topics listed above. At the time, 
only the LBL-TPC proposal was approved and became the STAR experiment. 
Eventually the MARS proposal evolved into the PHOBOS experiment (note 
that Phobos is a moon of the planet Mars) and PHENIX (should be spelled 
Phoenix) rose from the ashes of OASIS, Di-Muon, TALES and the other 
lepton focussed experiments. 

Eventually there were four approved experiments which have now been 
constructed and operated during the first year of RHIC running. BRAHMS, 
PHENIX, PHOBOS and STAR are briefly described below. These experi- 
ments have various approaches to study the deconfinement phase transition 
to the quark gluon plasma. The STAR experiment [7] concentrates on mea- 
surements of hadron production over a large solid angle in order to measure 
single- and multi-particle spectra and to study global observables on an event- 
by-event basis. The PHENIX experiment [8] focuses on measurements of 
lepton and photon production and has the capability of measuring hadrons 
in a limited range of azimuth and pseudo-rapidity. The two smaller experi- 
ments BRAHMS (a forward and mid-rapidity hadron spectrometer) [9] and 
PHOBOS (a compact multiparticle spectrometer) [10] focus on single- and 
multi-particle spectra. The collaborations, which have constructed these de- 
tector systems and which will exploit their physics capabilities, consist of 
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approximately 900 scientists from over 80 institutions internationally. In ad- 
dition to colliding heavy ion beams, RHIC will collide polarized protons to 
study the spin content of the proton [11]. STAR and PHENIX are actively 
involved in the spin physics program planned for RHIC. 

4.2. BRAHMS 

The BRoad Range Hadron Magnetic Spectrometer BRAHMS experiment 
is designed to measure and identify charged hadrons (tt^, K^, (p)) over a 
wide range of rapidity and transverse momentum for all beams and energies 
available at RHIC. Because the conditions and thus the detector requirements 
at mid-rapidity and forward angles are different, the experiment uses two 
movable spectrometers for the two regions. 




Figure 4. Layout of the BRAHMS detector. 

As shown in Figure 4, there is a mid-rapidity spectrometer to cover the 
pseudo-rapidity range < rj < 1.3 and a forward spectrometer to cover 
1.3 < rj < 4.0. The latter employs four dipole magnets, three time projec- 
tion chambers (TPC), and drift chambers. Particle identification is achieved 
with time-of-flight hodoscopes, a threshold Cherenkov counter, and one ring- 
imaging Cherenkov counter (RICH). The solid angle acceptance of the for- 
ward arm is 0.8 mstr. The mid-rapidity spectrometer has been designed for 
charged particle measurements for p < 5 GeV / c. The spectrometer has two 
TPCs for tracking, a magnet for momentum measurement, and a time-of- 
flight wall and segmented gas Cherenkov counter (GASC) for particle identi- 
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fication. It has a solid angle acceptance of 7 mstr. A set of beam counters and 
a silicon multiplicity array provide the experiment with trigger information 
and vertex determination. 

4.3. PHOBOS 




Figure 5. PHOBOS detector setup for the 2000 running period. 



The PHOBOS detector is designed to detect as many of the produced 
particles as possible and to allow a momentum measurement down to very 
low p±. The setup consists of two parts: a multiplicity detector covering 
almost the entire pseudo-rapidity range of the produced particles and a two 
arm spectrometer at mid-rapidity. Figure 5 shows the detector, including the 
spectrometer arms, the multiplicity and vertex array, and the lower half of the 
magnet. 

One aspect of the design is that all detectors are produced using a com- 
mon technology, namely as silicon pad or strip detectors. The multiplicity 
detector covers the range —5.4 < rj < 5.4, measuring total charged multi- 
plicity dN cn /d?7 over almost the entire phase space. For approximately 1% of 
the produced particles, information on momentum and particle identification 
will be provided by a two arm spectrometer located on either side of the 
interaction volume (only one arm was installed for the 2000 run). Each arm 
covers about 0.4 rad in azimuth and one unit of pseudo-rapidity in the range 
< T] < 2, depending on the interaction vertex, allowing the measurement 
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of p± down to 40 MeV/c. Both detectors are capable of handling the 600 Hz 
minimum bias rate expected for all collisions at the nominal luminosity. 

4.4. PHENIX 

The PHENIX experiment is specifically designed to measure electrons, muons, 
hadrons and photons. The experiment is capable of handling high event rates, 
up to ten times RHIC design luminosity, in order to sample rare signals such 
as the J/ip decaying into muons and electrons, high transverse momentum 
7r°'s, direct photons, and others. The detector consists of four spectrometer 
arms. Two central arms have a small angular coverage around central rapidity 
and consist of a silicon vertex detector, drift chamber, pixel pad chamber, ring 
imaging Cerenkov counter, a time-expansion chamber, time-of-flight and an 
electromagnetic calorimeter. These detectors allow for electron identification 
over a broad range of momenta in order to measure both low mass and high 
mass vector mesons. Two forward spectrometers are used for the detection 
of muons. They employ cathode strip chambers in a magnetic field and inter- 
leaved layers of Iarocci tubes and steel for muon identification and triggering. 
The overall layout of the PHENIX detector is shown in Figure 6. 



PHENIX Detector - Second Year Physics Run 




West Beam View East 

Figure 6. Shown is a beam view of the PHENIX two central spectrometer arms. Their is 
a axial field magnet in the middle. The detectors from the inner radius out are the multi- 
plicity and vertex detector (MVD), beam-beam counters (BBC), drift chambers (DC), pad 
chambers (PC 1-3), ring imaging cherenkov counter (RICH), time-expansion chamber (TEC), 
time-of-flight (TOF), and a Lead Glass and Lead Scintillator electro-magnetic calorimeter 
(PbSc, PbGl). 
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PHENIX Detector - Second Year Physics Run 




□ Installed □ Active 

Figure 7. Shown is a side view of the PHENIX detector including the two central spectrometer 
arms and the two muon spectrometers. The muons systems consist of cathode strip chamber 
muon trackers (MuTr) and muon identifiers (MuID) interleaved with layers of steel. 

One key feature of the PHENIX detector is the high rate capability of 
the data acquisition system (DAQ) and multi-level trigger architecture. These 
allow PHENIX to sample physics from RHIC collisions above the design 
luminosity of the machine. This high rate is crucial for studying rare leptonic, 
photonic and high p± processes. 

4.5. STAR 

The Solenoidal Tracker At RHIC (STAR) is a large acceptance detector capa- 
ble of tracking charged particles and measuring their momenta in the expected 
high multiplicity environment. It is also designed for the measurement and 
correlations of global observables on an event-by-event basis and the study 
of hard parton scattering processes. The layout of the STAR experiment is 
shown in Figure 4.5. The initial configuration of STAR in 2000 consists of 
a large time projection chamber (TPC) covering \rj\ < 2, a ring imaging 
Cherenkov detector covering \rj\ < 0.3 and A0 = 0.17T, and trigger detectors 
inside a solenoidal magnet with 0.25 T magnetic field. The solenoid provides 
a uniform magnetic field of maximum strength 0.5 T for tracking, momentum 
analysis and particle identification via ionization energy loss measurements 
in the TPC. Measurements in the TPC were carried out at mid-rapidity with 
full azimuthal coverage (Ac/> = 27r) and symmetry. A total of 1M minimum 
bias and 1M central events were recorded during the summer run 2000. 
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Figure 8. Schematic view of the STAR detector. 

Additional tracking detectors will be added for the run in 2001. These 
are a silicon vertex tracker (SVT) covering \rj\ < 1 and two Forward TPCs 
(FTPC) covering 2.5 < |r/| < 4. The electromagnetic calorimeter (EMC) will 
reach approximately 20% of its eventual — 1 < r] < 2 and A</> = 2ir coverage 
and will allow the measurement of high transverse momentum photons and 
particles. The endcap EMC will be constructed and installed over the next 2 
- 3 years. 

4.6. RHIC Spin Program 

The design of both the STAR and PHENIX experiments includes a polarized 
proton program to conduct studies of the spin structure of the proton. Crit- 
ical to this measurement is the identification of high transverse momentum 
photons and leptons. The STAR experiment is phasing in an electromagnetic 
calorimeter that will be crucial for such observations. In particular the RHIC 
experiments are ideally suited for measuring the gluon contribution of the 
proton spin. 
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5. Experimental Techniques 

5.1. Tracking Charged Particles (The STAR Example) 

The STAR experiment aims at the observation of hadronic observables and 
their correlations, global observables on an event-by-event base, and the mea- 
surement of hard scattering processes. The physics goals dictate the design 
of an experiment. Given the physics directions it is easy to summarize the 
general requirements: 

- Soft physics (100 MeV/c < p± < 1.5 GeV/c) 

• detection of as many charged particles as possible with high effi- 
ciency to provide high statistics for event-by-event observables and 
fluctuation studies 

• 2tt continuous azimuthal coverage for reliable event characteriza- 
tion 

• high tracking efficiency as close to the vertex as possible to contain 
the size of the experiment 

• adequate track length for tracking, momentum measurement and 
particle identification for a majority of particles 

• good two-track resolution providing a momentum difference reso- 
lution od a few MeV/c for HBT studies 

• accurate determination of secondary vertices for detecting strange 
particles (A, S, JT) 

— Hard physics (> 1.5 GeV/c and jets) 

• large uniform acceptance to maximize rates and minimize edge 
effects in jet reconstruction 

• accurate determination of the primary vertex in order to achieve 
high momentum resolution for primary particles 

• electromagnetic calorimetry combined with tracking and good mo- 
mentum resolution up to p± =12 GeV/ c to trigger on jets 

• segmentation of electromagnetic calorimeters which is consider- 
able finer than the typical jet size, i.e. jet radius r = \J dr\ 2 + d(j) 2 ~ 
1. 

The challenge is to find a detector concept which meets all these requirements 
with minimal costs. The detector of choice to solve the main tracking tasks 
was a large Time Projection Chamber (TPC) operated in a homogeneous 
magnetic field for continuous tracking, good momentum resolution and parti- 
cle identification (PID) for tracks below 1 GeV/c. The requirements not met 
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by the TPC needed to be covered by more specialized detectors such as a large 
acceptance electromagnetic calorimeter (hard processes), two Forward-TPCs 
(coverage of 77 > 1.7), and an inner silicon vertex tracker (better primary and 
secondary vertex measurement). For particle identification of high-p^ parti- 
cles detectors a small Ring Imaging Detector (RICH) and a Time-of-Flight 
(ToF) patch were added. Both detectors are not usable for event-by-event 
physics due to their small acceptance but allow to extend the PID capabilities 
for inclusive distributions. 

The need for an homogeneous field along the beam direction puts an 
stringent constraint on the design of the whole experiment. Only large so- 
lenoidal magnets are able to provide uniform fields of considerable strength 
(0.5 T). To keep down the costs the magnet cannot be too large which limits 
significantly the amount of "real estate" (i.e. detectors) it can house inside 
for tracking, PID, and calorimetry. The final magnet has coils with an inner 
radius of 2.32 m and a yoke radius of 2.87 m. The total length is 6.9 m. 
Note, that this concept is very different from the design of PHENIX where 
the axial-field magnet does not, or only weakly, constrain the dimensions of 
the required detectors. 



Figure 9. Peripheral Au+Au event recorded in the STAR TPC. Shown is the projection of all 
hits (points) and reconstructed tracks (solid lines) in the event onto the xy plane perpendicular 
to the beamline. 
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In the following we focus on the TPC, STARs main tracking detector, 
that essentially performs the role of a 3D camera with around 70 million 
pixel resolution. The TPC is divided into two longitudinal drift regions, each 
2.1 m long. Electrons created from track ionization drift in the longitudinal 
direction, along the TPC electric field lines, to the end-caps of the TPC. 
Each end-cap is instrumented with 70,000 pads. Each pad reads out 512 time 
samples. The position of the ionization charge in the readout plane provides 
the x and y coordinates of a space point along the particle trajectory while 
the arrival time of the charge allows to determine the original z position. The 
ionization pattern (chain of hits) of a traversing charged particle curved in 
the magnetic field allows the complete reconstruction of the particle trajec- 
tory and its 3-momentum. Fig. 9 shows the xy-projection (front view) of a 
low multiplicity event recorded in the STAR TPC. Each point represents one 
reconstructed hit, i.e. the local ionization charge created by one particle. The 
lines represent the reconstructed trajectories of the particles. Together with 
the timing information it is possible to also reconstruct the z-position of the 
hits and such the polar angle of the tracks as shown in Fig. 10. Sophisticated 
pattern recognition programs perform the task of reconstructing the particle 
trajectory using the measured positions of the measured hits. This procedure 
is commonly referred to as "tracking". In the following we discuss the pa- 
rameterization used to describe a track in the STAR TPC (and in any other 
detector with an homogeneous solenoidal field). 




Figure 10. Same event as shown in Fig. 9 but viewed from the side. The beamline runs 
horizontally from the left to the right. 
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5.1.1. Track Parameterization Momentum Determination 
The trajectory of a charged particle in a static uniform magnetic field with 
B = (0, 0, B z ) is a helix. In principle five parameters are needed to define 
a helix. From the various possible parameterizations we describe here the 
version which is most suited for the geometry of a collider experiment and 
therefore used in STAR. 

This parameterization describes the helix in Cartesian coordinates, where 
x, y and z are expressed as functions of the track length s. 



x(s) = xq H — [cos(<£o + h s k cos A) — cos <J>q] (1) 

K 

t/(s) = Vo + - [sin($ + hs k cos A) - sin <J? ] (2) 

K 

z(s) = zq + s sin A (3) 

where: s is the path length along the helix 
x o j yo j zo is the starting point at s = sq = 
A is the dip angle 
n is the curvature, i.e. k = 1/R 

B is the z component of the homogeneous magnetic field (B = (0, 0, B z )) 
q is charge of the particle in units of positron charge 

h is the sense of rotation of the projected helix in the xy-plane, i.e. h = 
— sign(qB) = ±1 

$o is the azimuth angle of the starting point (in cylindrical coordinates) with 
respect to the helix axis ($o = \£ — hir/2) 

* is the arctan(dy/dx) s= o, i.e. the azimuthal angle of the track direction at 
the starting point. 

The meaning of the different parameters is visualized in Fig. 11. 

The circle fit in the xy-plane gives the center of the fitted circle (x c , y c ) 
and the curvature k = 1/R while the linear fit gives zq and tan A. The phase 
of the helix (see Fig. 1 1) is defined as follows: 

<D = arctan ( ^LZ^") (4) 



,20 - X c , 

The reference point (xq, yo) is then calculated as follows: 

cos $o 

xo = x c H (5) 

K 

yo = 2/cH (6) 

K 

and the helix parameters can be evaluated as: 

* = $ Q + hir/2 (7) 



corsica.tex; 8/02/2008; 21:06; p. 20 



Experimentalists and Theorists 



21 




s>0 



y ii 



t 



x 



s 



xy 



Figure 11. Helix parameterization: shown on the left is the projection of a helix on the xy 
plane. The crosses mark possible data points. The right plot depicts the projection of a helix 
on the sz plane. For the meaning of the various parameters see text. 



where k is the curvature in [m ], B the value of the magnetic field in [Tesla], 
c the speed of light in [m/ns] 0.3) and p± and p z are the transverse and 
longitudinal momentum in [GeV/c]. 

Once the track momenta and parameters are determined, one can use 
the above parameterization to extrapolate the TPC tracks to other detectors, 
e.g. outwards to the electromagnetic calorimeter or RICH detector, or inwards 
to the silicon vertex tracker (SVT). Once the referring hits are found they are 
added to the track and the parameters are re-evaluated. The charges of all hits 
along a track can be used to calculate its energy loss (dE / dx) and together 
with its known momentum can be used to determine the particle mass. 

5.2. Electron Identification (the PHENIX example) 

As we have seen in the physics motivation section, there are many reasons 
to measure electrons over a broad range in transverse momentum. In the 
original Letters of Intent there were many different proposals for detector 
technologies for measuring electrons, and many of these capabilities were 
combined into the PHENIX experiment. The PHENIX experiment has two 
central spectrometer arms each covering 90 degrees in (j> and |^| < 0.35. 

Electrons and also charged hadrons that are produced near mid-rapidity 
are bent in an axial magnetic field over a radial distance of approximately 



Pz 



c q B I k 
pj_ tan A 



(8) 
(9) 




(10) 
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Figure 12. A schematic diagram of the PHENIX central spectrometer magnetic reconstruction 
technique is shown. 



two meters after which the aperture is relatively field free. There is a large 
aperture drift chamber that measures the projective trajectory of the charged 
particle tracks in the field free region with multiple wire layers oriented in 
the x direction (giving maximum resolution in the bend plane of the field) 
and also in u and v direction (stereo planes for pattern recognition). The 
drift chamber is augmented by a series of moderate resolution pad chambers, 
which yield three dimensional space point along the particle's track and aid 
significantly in pattern recognition and non-vertex background rejection. As 
shown in Figure 12 the particle track is characterized by an angle (a) in the 
bend plane of the magnetic field. From this vector and the assumption that 
the track originated at the Au-Au collision vertex, the rigidity of the track is 
determined. The rigidity is the momentum divided by the particle charge, and 
with the assumption of a Z = ±1 particle, the momentum is known. 

Typically the momentum resolution has two components as shown be- 
low. 

Sp/p = (« 1%) + (« 1%) x p[GeV/c] (1 1) 

The first term, typically of order 1% is due to the multiple scattering of the 
charge particle in material before and in the tracking devices. The second 
term which scales with the particle's momentum is related to the finite spatial 
resolution of the detector. The PHENIX detector has been designed with a 
minimum of inner region material including a Beryllium beam pipe, which 
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has a low Z value to minimize multiple scattering while maintaining integrity 
for the vacuum. All four experiments at RHIC have Be beam-pipes in the in- 
teraction region. The PHENIX design resolution from the drift chamber is of 
order 0.5% at p± « 0.2 GeV and increases linearly with p± above 0.7 GeV. 
There is a great deal of work involved in calibrations and wire alignement to 
achieve these resolution values. 

It is notable that PHENIX uses an axial field magnet for the central spec- 
trometers as opposed to a solenoidal field. There are substantial advantages 
and disadvantages to this choice. One major consideration is that if the entire 
experiment is contained inside a large solenoid, as in the STAR configuration, 
the magnetic field is present throughout the volume. This is advantageous 
for the large volume Time Projection Chamber (TPC) of STAR where they 
can observe the curvature of the charged particle tracks. However, the cost 
of these magnets increases steeply as one increased the desired outer radius. 
PHENIX for the purposes of hadron particle identification wanted to have a 
Time-of-Flight scintillator wall approximately 5 meters from the interaction 
vertex and for electron and photon identification a good energy resolution and 
high granularity electromagnetic calorimeter behind that. In keeping within 
the experimental budget it would not have been possible to have a solenoidal 
magnet with an outer radius greater than 5.5 meters. The ALICE experiment 
being built for the CERN-LHC program is fortunate to be able to re-use the 
very large solenoid from the L3 experiment at LEP. One disadvantage of the 
PHENIX choice is that the pole tips of the axial field magnet are close to the 
interaction point and create substantial "shine", particles scattering off the 
poles into the detector aperture. The overall choice for the field configuration 
is an important starting point for many detector designs. 

Now that we have characterized the particle's momentum vector, we 
must discriminate all of the charged pions from our interesting electrons. 
The first detector in PHENIX that is employed is a Ring Imaging Cherenkov 
detector (RICH). It is a large gas volume detector with a thin mirror plane 
for reflecting Cherenkov light onto an array of photo-multiplier tubes (PMT) 
that are situated off to the side of the spectrometer acceptance. The radiator 
gas used is either ethane with an index of refraction n= 1.00082 or methane 
with n= 1.00044. Requiring more than three PMT hits yields almost 100% 
efficiency for electrons and rejects pions with p± 4 GeV at the level 
of < 10 3 in a single track environment. The detector, requiring only three 
PMT's, does not reconstruct a ring radius for further particle characterization, 
but rather is used as a threshold detector only. 

In a multiple track collisions (remember of order 5000 charged particles 
are being produced) one can incorrectly match a track to a RICH signal. 
Further electron identification is provided by the electromagnetic calorimeter 
(EMCal). The calorimeter is composed of both Lead Scintillator (PbSc) and 
Lead Glass (PbGl) modules, with the later being originally used in the WA98 
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experiment at the CERN-SPS heavy ion program. The calorimeter has a ra- 
diation length of I8X0 and m 16Xq for the PbSc and PbGl respectively, 
but does not fully interact and contain hadronic showers. Thus, an electron 
or photon incident on the calorimeter deposits most of its energy, while a 
hadron (eg. charged pion) has a large probability to pass through the module 
depositing a small dE/dx minimum ionizing radiation energy. Even when 
the pion suffers an inelastic collision in the calorimeter, only a fraction of 
its energy is contained and measured in the PMT at the back of the module. 
The excellent energy resolution m 5 — 8Vo/y/E(GeV) and high granularity 
give precise electron identification and pion rejection by requiring the energy 
match the measured momentum (E/p « 1). However, the background rejec- 
tion degrades when the particle momentum is low and hadronic showers have 
a higher probability to match the measured momentum. 

Therefore, an additional detector is necessary to help identify low p± 
electrons for the crucial low mass vector meson physics. For this purpose, 
PHENIX uses a Time Expansion Chamber (TEC) that samples energy loss in 
a gas radiator. The TEC determines the particle species using dE/dx infor- 
mation. It has a rejection of e/ir « 5% for particles with a momentum p = 500 
MeV/c with P10 gas and 2% with Xe gas, which is much more expensive. 
It is the combination of all these detectors that allow for efficient electron 
measurements with a minimum of pion contamination. Thus, the challenge 
of electron identification over a broad range in p± is met. 

5.3. Muon Identification (the PHENIX example) 

Direct muons and hadrons decaying into muons require a rather different 
experimental approach to measuring electrons, photons and other hadrons. 
Muons interact with a low cross section in material, and are easiest to identify 
by placing steel or other material in the particle path and removing all other 
particles through interaction. Then a detector placed after the steel should 
measure a clean muon sample. 

However, the muons must have a large enough energy to penetrate the 
steel without stopping due to ionization energy loss dE/dx. In the central 
rapidity region at a collider, the muons from low p± J/tp decays make it 
impossible to have enough steel to range out other hadrons effectively, while 
allowing the muons to pass through. In particle experiments that focus on 
high p± muons, the detector can consist of a similar one to that described 
above for electrons. After the calorimeter, one can have some steel absorber 
(often the return in the magnet steel) and then have a muon identifier detector 
outside of that. This design does not work well for heavy ion physics. First, 
it restricts the measurable p± range at too high a value and also has a large 
background from low p± pions decaying into muons. 
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Another option is to measure muons are forward rapidity where they 
have a substantial momentum in the longitudinal direction. PHENIX mea- 
sures muons in the forward and backward pseudo-rapidity regions (at angles 
of 10-35 degrees from the beam line). The detector consists of a brass and 
steel absorber to range out hadrons followed by a cathode strip chamber muon 
tracking device that measures the particle's bend in a magnetic field. There is 
a delicate balance in the amount and type of absorber material used. Too much 
material and the multiple scattering reduces the resolution which is important 
for cleanly separating states such as J/tp and ij/, and too little material in 
which case the particle occupancy in the tracking device is too large. After 
the muon tracking device there is more steel absorber with Iarocci tube muon 
identifiers interspersed. The coverage is from 1.1 < < 2.4 and the muons 
must have E > 2.1 GeV in order to penetrate the absorber material. The 
identifier detectors also provide the necessary trigger information to sample 
the muons from the high luminosity RHIC collisions. 

The most substantial background in measuring J/ ip and D mesons are 
muons from pions and kaons that decay before they hit the brass nose-cone 
absorbers. There is a competing requirement in PHENIX in that one wants 
the absorber as close to the interaction vertex as possible to reduce this decay 
contribution, but they need to be far enough apart for there to be a good accep- 
tance for electrons, hadrons and photons in the two central arm spectrometers. 
There are many benefits to the comprehensive design of PHENIX, but there 
are definite drawbacks as well. Of the two muon arms, the south muon arm 
was completed for running in Run II at RHIC, and the north arm will be 
complete for Run III. 



6. Physics Results 

In this chapter we present the physics results from the Run I data and a 
preview of expected results from Run II. 

6.1. Global Observables 

The first result with a measurement from all four RHIC experiments is the 
charged particle multiplicity [12, 13, 14, 15]. The four experiments' results 
for central (small impact parameter) collisions are in excellent agreement and 
are shown in Figure 13. The multiplicity rises more sharply as a function of 
center-of-mass energy in heavy ion collisions than in p+p and p+p collisions, 
which is attributed to the increased probability for hard partem scattering in 
the thick nuclear target seen by each partem. 

We expect the charge particle yield to increase for collisions of larger 
nuclei. However, at low x values, the high density of gluons may in fact 
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Figure 13. Charged particle multiplicity measurement from all four RHIC experiments is 
shown for Au-Au collisions at ^/s NN = 130 GeV. Also shown are data for p + p and p + p 
collisions. A model of heavy ion collisions HIJING is shown for comparison. 

saturate due to gluon fusion processes. The contribution to the yield from 
hard processes should exhibit point-like scaling (scaling with the number of 
binary collisions) and would thus scale as A 4 / 3 . However, parton saturation 
depends upon the nuclear size and would limit the growth of the number of 
produced partons as A 1 / 3 . If present, this initial parton saturation would limit 
the hard process contribution to the total charged particle multiplicity. 

Only one nuclear species (Au) was accelerated in Run I at RHIC. Thus, 
rather that changing the mass number A directly we control the collision 
volume by varying the centrality or the number of participating nucleons 
for Au-Au collisions. Shown in Figure 14 are the published results from 
the PHENIX [13] and PHOBOS [16] experiments for the number of charged 
particles per participant nucleon pair as a function of the number of partic- 
ipating nucleons. The number of participating nucleons is determined in a 
slightly different manner by the different experiments. However, the general 
method is to calibrate the number of spectator nucleons (= 2 x A - participant 
nucleons) using a measurement of spectator neutrons in a set of zero degree 
calorimeters that are common to all experiments. By correlating the number 
of forward neutrons to the number of charged particle produced in the large 
pseudo-rapidity region, the event geometry can be understood. 

In Figure 14 one can also see theory comparisons that indicate that a 
model including parton saturation (EKRT [17]) fails to agree with the more 
peripheral data. Results from the HIJING model [18] are also shown which 
does not include parton saturation and thus has a more continuous rise in 
the particle multiplicity. Since saturation phenomena are only likely to have 
observable consequences for large collision volumes, it is not possible with 



corsica.tex; 8/02/2008; 21:06; p. 26 



Rxnerimenralists and Theorists 



27 




Z i I , I , , i 

100 200 300 400 

N part 

Figure 14. PHENIX and PHOBOS results for dN ch /dr)\ v=0 / ±N part as a function of N part . 
The hashed and solid bands indicate the systematic errors for the two experimental results. The 
data point for pp with two participants is shown for comparison. Also theoretical predictions 
from the HIJING and EKRT models are shown. 



present systematics to rule out the saturation picture for the most central 
collisions. 

In order to better test the saturation picture lighter ion, smaller A, colli- 
sions will be studied in Run II. In addition, heavy flavor (charm and bottom) 
and Drell-Yan production should be a sensitive probe to the initial parton 
density. Another proposal is that by varying the collision energy and keeping 
the nuclear geometry the same one can get a better handle on systematics and 
test scenarios dependent on the coupling constant and the saturation scale. 
The physics of parton saturation and color glass condensates is at the fore- 
front of theoretical development in the field. Many recent developments were 
discussed in this workshop and will be described in other contributions. 

In addition to the initial parton density, the energy density is of great 
interest. There are published results estimating the initial thermalized energy 
density achieved in these collisions. Bjorken originally derived a formula, 
shown in Eqn. 12, relating the measured transverse energy per unit rapidity 
to the thermal energy density [19]. 

1 1 dE ± 

€b j = ~m T~ ( lT> 

J ttK z ct ay 

It should be noted that there is a trivial factor of two error in the original 
reference that is corrected here. This formulation assumes a boost invari- 
ant expanding cylinder of dense nuclear matter and a thermalization time r. 
There are two important assumptions in this particular formulation. The first 
is the boost invariant nature of the collision. There are recent preliminary 
measurements from STAR and PHOBOS that indicate the distribution of 
particles is relatively flat over ±2 units of pseudo-rapidity. However, shown 
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Figure 15. Plotted is the n /tt + and p/p ratio as a function of rapidity from the BRAHMS 
experiment for central collisions. 
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Figure 16. The PHENIX experiment result is shown for dE±/dri/(0.5Np) at r\ — as a 
function of the number of participating nucleons. Also shown in the result from experiment 
WA98 at the lower energy CERN-SPS. 



in Figure 15 is the measured distribution of p/p from the BRAHMS experi- 
ment [20]. This indicates the the system is already changing at y « 2, though 
it is not clear that this is enough to invalidate the energy density formulation. 
The second question is what is the relevant thermalization time r. 
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The PHENIX experiment has published [21] the transverse energy dis- 
tribution for minimum bias Au-Au collisions. For the 5% most central events, 
the extracted transverse energy < dEj_ /drj > \ v= o = 503 ± 2 GeV. Shown in 
Figure 16 is dE±/dr]/(0.5N p ) versus the number of participating nucleons. 
One sees a similar increase in transverse energy as was seen in the charged 
particle multiplicity yield. 

The canonical thermalization time used in most calculations is r = 1 
fm/c, that yields an energy density of 4.6 GeV/fm 3 , which is is 60% larger 
than measured at the CERN-SPS. In addition, it is believed that the density is 
substantially higher due to the potentially much shorter thermalization time 
in the higher partem density environment. If one achieves gluon saturation the 
formation time is of order 0.2 fm/c and gives an estimated energy density 
of 23.0 GeV/fm 3 . There are even estimates of over 50 GeV/fm 3 , but they 
assume a very large drop in the final measured transverse energy due to work 
done in the longitudinal expansion of the system. All of these estimates are 
above the energy density of order 0.6 — 1.8 GeV/fm 3 corresponding to the 
phase transition temperature 150 — 200 MeV. 

6.2. Elliptic Flow 

The azimuthal anisotropy of the transverse momentum distribution for non- 
central collisions is thought to be sensitive to the early evolution of the sys- 
tem. The second Fourier coefficient of this anisotropy, V2, is called elliptic 
flow [22]. It is an important observable since it is sensitive to the re-scattering 
of the constituents in the created hot and dense matter. This re-scattering 
converts the initial spatial anisotropy, due to the almond shape of the over- 
lap region of non-central collisions, into momentum anisotropy. The spatial 
anisotropy is largest early in the evolution of the collision, but as the sys- 
tem expands and becomes more spherical, this driving force quenches itself. 
Therefore, the magnitude of the observed elliptic flow reflects the extent 
of the re-scattering at early time [23]. The time evolution of the transverse 
energy density profile is schematically depicted in Figure 17 where the solid 
lines represent surfaces of constant energy density. The pressure in the system 
is highest in direction of the reaction plane (largest energy density gradient) 
which causes the elliptic anisotropy. 

Elliptic flow in ultra-relativistic nuclear collisions was discussed as early 
as 1992 [24] and has been studied intensively in recent years at AGS [25, 26], 
SPS [27, 28, 29] and now at RHIC [30] energies. The studies at the top AGS 
energy and at SPS energies have found that elliptic flow at these energies is 
in the plane defined by the beam direction and the impact parameter, V2 > 0, 
as expected from most models. 

The STAR detector is especially suited to study elliptic flow due to its 
azimuthal symmetry, large coverage, and its capability of tracking charged 
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Figure 17. Schematic view of a evolution of the transverse en ergy density profile (indicated 
by constant energy density con tours spaced by) and of the flow velocity field (indicated by 
arrows) for Pb+Pb collisions at impact parameter b=7.0 fm. The four panels show snapshots 
at times r — To = 3.2, 4.0, 5.6, and 8.0 fm/c. At these times the maximal energy densities in 
the center are 5.63, 3.62, 1.31 and 0.21 GeV/fm 3 , respectively. The figure is taken from [32]. 

particles down to very low p±. Even in peripheral events there are sufficient 
tracks available to divide the event in two 'subevents' of which one is used 
to measure, or better estimate, the event plane and the other to correlate the 
particles in it in order to derive V2- While BRAHMS is not able to study 
elliptic flow because of its small acceptance, PHOBOS has some capabilities 
to measure vi but only integrated over all p±. Because of its restricted az- 
imuthal coverage the PHENIX collaboration follows a different approach to 
reconstruct v-i by studying the A<p correlation between particles thus circum- 
venting the event plane determination. However, work is still in progress and 
in the following we concentrate on the more direct event plane method used 
in STAR. 

The flow analysis method involves the calculation of the event plane 
angle, which is an experimental estimator of the real reaction-plane angle. 
The second harmonic event plane angle, ^2, is calculated for two sub-events, 
which are independent subsets of all tracks in each event. Figure 18 shows 
the results for the correlation between the sub-events for the first and second 
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Figure 18. The correlation between the event plane angles determined for two independent 
sub-events. The correlation is calculated for the first harmonic (n=l) and the second harmonic 
(n=2). 



harmonic as a function of centrality [30]. The peaked shape of the centrality 
dependence of (cos[2(^ a — is a signature of anisotropic flow. However, 
the correlation between the sub-events may not be due entirely to anisotropic 
flow. To estimate the magnitude of non-flow effects one can use the sub- 
events in three different ways: 

1. Assigning particles with pseudo-rapidity 7] < + e to one sub-event 
and particles with i] > + e to the other. Short range correlations, such 
as Bose-Einstein or Coulomb, are to a large extent eliminated by the 2e 
"gap" between the two sub-events. 

2. Dividing randomly all particles into two sub-events, sensitive to all non- 
flow effects. 

3. Assigning positive particles to one sub-event and negative particles to the 
other, allowing an estimation of the contribution from resonance decays. 

Studies have shown that the results from all three methods are for the 
central and mid-peripheral events very similar. For the most peripheral events 
the results vary among the methods by about 0.005. Not all non-flow contri- 
butions might be known and the effects of others, such as jets, are difficult 
to estimate because of their long-range correlation. In order to estimate the 
systematic uncertainty due to the effects of jets, one can assume that jets con- 
tribute at the same level to both the first and second order correlations. Taking 
the maximum observed positive first order correlation, as being completely 
due to non-flow will reduce the calculated V2 values. 
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Figure 19. The integrated elliptic flow signal, V2, with the estimated systematic uncertainties 
as measured by STAR. 
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Figure 20. Excitation function of V2 from top AGS to RHIC energies. 



Figure 19 shows the final V2 integrated over all p±_ as a function of 
centrality. The statistical uncertainties are smaller than the markers and the 
uncertainties shown are the systematic uncertainties due to this estimated 
non-flow effect. 

Figure 20 shows the maximum V2 value as a function of collision energy. 
It rises monotonically from about 0.02 at the top AGS energy [25], 0.035 at 
the SPS [28] to about 0.06 at RHIC energies [30]. This increasing magni- 
tude of the integrated elliptic flow indicates that the degree of thermalization, 
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which is associated with the amount of re-scattering, is higher at the higher 
beam energies. However, interpretation of the excitation function has to be 
done with care. The V2 values used here are the maximum values as a function 
of centrality for each energy. The centrality where V2 peaks can change as a 
function of beam energy, indicating different physics [31]. 
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Figure 21. V2 as a function of p±, as measured by STAR in v /s^7 = 130 GeV Au-Au 
collisions. 



The differential anisotropic flow is a function of r/ and p±. Figure 21 
shows V2 for charged particles as a function of p± for a minimum bias event 
sample. Mathematically, the v 2 value at p±=0, as well as its first derivative, 
must be zero, but it is interesting that v 2 appears to rise almost linearly with 
p± starting from relatively low values of pj_This is consistent with a stronger 
"in-plane" hydrodynamic expansion of the system than the average radial 
expansion. 

Comparing to estimates based on transport cascade models, one finds 
that elliptic flow is under predicted by a factor of more than 2. Hydrodynamic 
calculations [32] for RHIC energies over predict elliptic flow by about 20- 
50%. This is just the reverse of the situation at the SPS where cascade models 
gave a reasonable description of the data and hydrodynamic calculations were 
more than a factor of two too high. Also in contrast to lower collision en- 
ergies, the observed shape of the centrality dependence of the elliptic flow 
is similar to hydrodynamic calculations and thus consistent with significant 
thermalization which is one of the most striking results from the initial round 
of RHIC results. 
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6.3. Two-Particle Interferometry (HBT) 



The study of small relative momentum correlations, a technique also known 
as HBT [33] interferometry, is one of the most powerful tools to study com- 
plicated space-time dynamics of heavy ion collisions [34]. It provides cru- 
cial information which helps to improve our understanding of the reaction 
mechanisms and to constrain theoretical models of the heavy ion collisions. 
Interpretation of the extracted HBT parameters in terms of source sizes and 
lifetime is more or less straightforward for the case of chaotic static sources. 
In the case of expanding sources with strong space-momentum correlations 
(due to flow, etc.) the situation is more difficult, but the concept of length of 
homogeneity [35] provides a useful framework for the interpretation of data. 

The dependence of the pion-emitting source parameters on the trans- 
verse momentum of the particle pairs (Kt) and on centrality can in principle 
be measured by all RHIC experiments with high statistics. For more detailed 
analysis as for example event-by-event HBT, HBT radii versus reaction plane, 
and the correlation of HBT results with other observables can only be per- 
formed by STAR due to its large acceptance and azimuthal coverage. These 
studies, however, are still in progress and it is by far too early to discuss 
them here. In the following we show results from the STAR experiment [36] 
that performed a multi-dimensional analysis using the standard Pratt-Bertsch 
decomposition [37] into outward, sideward, and longitudinal momentum dif- 
ferences and radius parameters. The data are analyzed in the longitudinally 
co-moving source frame, in which the total longitudinal momentum of the 
pair (collinear with the colliding beams) is zero. 

As expected, larger sizes of the pion-emitting source are found for the 
more central {i.e. decreasing impact parameter) events, which in turn have 
higher pion multiplicities. This source size is observed to decrease with in- 
creasing transverse momentum of the pion pair. This dependence is similar to 
what has been observed at lower energies and is understood to be an effect of 
collective transverse flow. Shown in Figure 22 is the coherence parameter A 
and the radius parameters R ou t, R s jde> and Riong obtained in the analysis. Also 
shown are values of these parameters extracted from similar analyses at lower 
energies. All analyses are for low transverse momentum (~ 170 MeV/c) neg- 
ative pion pairs at mid-rapidity for central collisions of Au + Au or Pb + Pb. 
From Figure 22 the values of A, R ou t, R s ide, and Rzong extend smoothly from 
the dependence at lower energies and do not reflect significant changes in 
the source from those observed at the CERN SPS energy. One of the biggest 
surprises is that the anomalously large source sizes or source lifetimes pre- 
dicted for a long-lived mixed phase [38] have not been observed in this study. 
Preliminary results of the HBT analysis by the PHENIX Collaboration [39] 
agree with the STAR results within error bars. 
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Figure 22. Compilation of results on two-particle correlation (HBT) parameters from mea- 
surements using central collisions of Au + Au at the BNL-AGS, Pb + Pb at the CERN-SPS 
and Au + Au data from the STAR experiment at RHIC. Plotted are the coherence parameter 
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One of the big puzzles, however, is the magnitude and the tranverse 
momentum (Kt) dependence of the ratio of R ou tlR s ide which contradicts 
all model predictions [38, 40]. These model calculations predict the ratio to 
be greater than unity due to system lifetime effects which cause Rout to be 
larger than R s id e - They also predict that the ratio increases with Kt- Such an 
increase seems to be a generic feature of the models based on the Bjorken- 
type, boost-invariant expansion scenario. Hence, it was surprising to see that 
the experimentally observed ratio is less than unity and is decreasing as a 
function of Kt- Currently, it is far from clear what kind of scenario can lead 
to such a puzzling Kt dependence. 

6.4. Hard Process Probes of the Plasma 
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Figure 23. Preliminary PHENIX invariant multiplicity of identified 7r as a function of trans- 
verse momentum are shown for peripheral and central collisions. Comparison with theoretical 
calculations with and without parton energy loss are also shown. 



Jet processes and their associated hadronic fragmentation provide one 
of the most exciting probes of the color deconfined plasma. The PHENIX 
experiment has measured the distribution of identified 7r° for both central 
and peripheral Au-Au collisions as shown in Figure 23 [41]. The peripheral 
results appear to be in good agreement within systematic errors of an ex- 
trapolation from pp collisions scaled up by the number of binary collisions 
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Figure 24. Invariant multiplicity of charged hadrons as a function of transverse momentum 
(top). Ratio of unidentified charge hadrons per calculated binary collision from Au + Au cen- 
tral collisions to those from p + p(p) collisions extrapolated to ^/s NN =130 GeV as a function 
of transverse momentum (bottom). 



expected in this centrality class. However, the central collision results show 
a significant suppression in the 7r° yield relative to this point like scaling 
expected for large momentum transfer parton-parton interactions. If the cre- 
ated fireball in RHIC collisions is transparent to quark jets, then we expect 
the yield of high pj_ hadrons to obey point-like scaling and equal the pp (or 
equivalently pp) distribution scaled up by the number of binary NN colli- 
sions, or equivalently by the nuclear thickness function Taa- This is not what 
is observed. A more sophisticated calculation [42] yields the same qualitative 
conclusion. 

The STAR experiment has recently submitted for publication [43] the 
p± spectra for unidentified negatively charged hadrons in central Au + Au 
collisions as shown in Figure 24. Also shown are the equivalent spectra from 
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Figure 25. PHENIX preliminary results for unidentified charged hadron invariant multiplicity 
as a function of transverse momentum. 

experiment NA49 at the CERN-SPS at yfs = 17 GeV and from UA1 in pp 
at y / s^" = 200 GeV. The STAR spectra is then divided by the spectra from 
pp scaled by the number of binary collisions, and the result is shown in the 
lower panel of Fig 24. At low transverse momentum the particle production 
is dominated by soft interactions which scale with the number of wounded 
nucleons as indicated by the line at 0.2. The rise from 0.2 as a function of 
p± certainly has a large contribution from hydrodynamic flow that will push 
particles to higher transverse momentum in central Au-Au collisions. 

The PHENIX experiment has shown preliminary results extending out 
further in transverse momentum. Preliminary results for six centrality classes 
are shown from PHENIX in Fig 25. STAR also has preliminary results for 
central collisions extending out to p± > 5 GeV that are in reasonable agree- 
ment with the PHENIX results. If one takes the ratio of the central spectra 
to the unidentified spectra in pp collisions scaled up by the number of binary 
collisions one gets a ratio Raa as shown in Fig 26. It needs to be noted 
that there is no pp data at -Js^ = 130 GeV and thus an extrapolation to 
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Figure 26. Ratio of unidentified charge hadrons per calculated binary collision from Au + Au 
central collisions to those from p + p(J5) collisions (extrapolated to y/s^ = 130 GeV) as a 
function of transverse momentum (GeV). The solid line is the systematic error band on the 
ratio. The dashed line is the average result from experiments at the lower energy CERN-SPS. 



that energy is done to calculate Raa- This extrapolation is included in the 
systematic error band, and should be reduced when both experiments measure 
the spectra in pp in Run II. 

There are many important physics points to understand in these results. 
The ratio appears to stay below one, although that is a marginal conclusion 
with the present systematic errors. However, this is certainly in qualitative 
agreement with a parton energy loss scenario, as also seen in the observed 
suppression in the PHENIX 7r° spectra. In contrast, the CERN-SPS results 
show an enhancement that has been attributed to the Cronin effect, or initial 
state parton scattering that gives a kx kick to the final transverse momentum 
distribution. This expected enhancement makes the suppression seen at RHIC 
all the more striking. 

There are a number of open questions that must be considered before 
drawing any conclusions. The most basic is that these p± values are low rel- 
ative to where one might have confidence in the applicability of perturbative 
QCD calculations. In addition, the separation between soft and hard scale 
physics is blurred in this p± range, and in fact the CERN-SPS ratio Raa has 
also been explained in terms of hydrodynamic boosting of the soft physics to 
higher p±. The preliminary results from PHENIX on the ratios of 7r/K/p(p) 
in the middle of this p± range look more like soft physics than a parton frag- 
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mentation function in vacuum. One additional point of concern is that these 
models of energy loss assume that the partem exits the collision region before 
finally fragmenting into a jet of forward hadrons. Thus the final hadronization 
takes place in vacuum. In the p± range of these early measurements, that 
conclusion is not so clear. The partem is traveling through the medium with 
various kr scatters, and if it hadronizes inside a bath of other particles, the 
leading hadrons may be slowed down by inelastic collisions with co-moving 
pions. Lastly, the point-like scaling is known to be violated due to the nuclear 
shadowing of parton distribution functions. These nuclear modifications are 
known to reduce the pdf for quarks of order 20% for x « 10~ 2 ; however, 
the shadowing for gluons is not currently measured. The calculations of [42] 
have included modeling of this shadowing, but must be viewed with caution 
at this time. These points need further theoretical investigation. In addition, 
as will be discussed in the next section, many of these concerns are reduced 
when the measurements extend to much higher transverse momentum. 



7. Future Measurements 

In the following chapter, we discuss a few select topics within our areas 
interest that have exciting results expected in the near future. 

7.1. Charm and Bottom 

The measurement of open charm and bottom in relativistic heavy ion collions 
is both an extreme experimental challenge and rich with physics informa- 
tion. First, the measurement of quarkonium states such as the J/ip require 
a comparison measure of the original cc production to determine the effect 
of color screening. Also, the total charm production is sensitive to the initial 
gluon density in the incoming nuclei and is thus sensitive to any shadowing of 
the gluon distribution function and may even comment on the possible color 
glass condensate postulated to describe the phase space saturated gluon distri- 
butions in the highly Lorentz contracted nuclei. Lastly, recent predictions of 
charm quark energy loss in tranversing a hot partonic medium have generated 
much interest. Now for the difficult part. The best was to measure charm via 
D mesons is either via direct reconstruction from its ir + K decay mode or 
via the semi-leptonic decay K + e + v e . The combinatoric background in 
the purely hadronic channel are close to overwhelming and the semi-leptonic 
decay cannot be completely reconstructed. In particle physics experiments 
measuring the decay products with a few micron displaced vertex from the 
collision vertex allows for a dramatic reduction in the combinatoric back- 
ground. However, the level of silicon detector technology was not advanced 
enough for this very high multiplicity environment at the time of the RHIC 
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Figure 27. Simulation of the transverse momentum spectrum for mid-rapidity single electrons 
and positrons. The top curve is the sum of all contributions, the two lower curves are from 
Dalitz decays and photon conversions. The subsequent three grey shaded areas are from charm 
D mesons, beauty B mesons, and Drell-Yan. 

detector designs. This is now being discussed as a possible upgrade to the 
experiments. 

One promising way of determining the charm production is through the 
measurement of single electrons. Shown in Figure 27 is a simulation of sin- 
gle electrons as a function of transverse momentum. The top curve in the 
sum of all contributions. The next two curves that dominate the spectra at 
low p± < 1.0 GeV and from the Dalitz decays of pions and rj and from 
photon conversions. The next two contributions are from charm D meson and 
beauty B meson decays. Charm yields « 50% of the counts at p± = GeV, and 
beauty yields > 50% of the counts above p± > 3.5 GeV. The lowest curve is 
the contribution from Drell-Yan which never has a major contribution to the 
single electrons. 

PHENIX has made a preliminary measurement of the single electron 
transverse momentum spectra from the limited statistics in Run I as shown 
in Figure 28. The analysis of these results is proceeding and implications on 
charm production are forthcoming. 
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Figure 28. PHENIX preliminary results for unidentified charged hadron invariant multiplicity 
as a function of transverse momentum. 



Additional handles on heavy flavor production can be had with the mea- 
surement of correlated leptons. For example, electron-muon pairs at large 
relative momentum (Q 2 or Mi nv ) have a substantial contribution from cc and 
bb pairs. For example, the cc can fragment into DD followed by the decays 
D — ► K + e + v e and D — > K + \i + v^. Although the initial Q 2 of 
the cc pair is significantly modified when measured as a Q 2 of the e/x pair, 
there is enough information to attempt to extract a total charm cross section 
and maybe something about the initial Q 2 of the cc. This measurement has the 
advantage over e + e~ and /i, + n~ pairs in that e/x pairs are free from Drell-Yan 
and thermal contributions. 

One of us (J.N.) has used PYTHIA 6.0 to estimate the rate of e/x pairs 
into the PHENIX acceptance (Central arms for the electron and South muon 
arm only for the //). PYTHIA has been run with a charm quark mass of 1.5 
GeV/c 2 and < {k T ) = 1.5 GeV/c. _ 

We show in Figure 29 the distribution of Q 2 for all cc, DD, e\i and e/i 
pairs accepted by PHENIX. One can see that the charm mesons carry most 
of the information from the cc pair. For the e/i pair the correlation with the cc 
pair Q 2 is substantially washed out and of much lower slope since the kaon in 
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Figure 29. Q 2 in GeV 2 for cc (black line), DD (red line), efi (blue line), and efi accepted by 
PHENIX (red fill). The vertical scale is arbitrary between the different curves. 
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the D decay takes away a large fraction of the original charm momentum that 
is not measured. It should be noted that the modeling of the fragmentation of 
the charm quark will be a source of systematic error, whereas the blurring of 
the Q 2 from the efi decay kinematics can be modeled exactly. 

We show the invariant mass distribution of the e/x pairs into the PHENIX 
acceptance in Figure 30. Also shows are three mass distribution with electron 
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energy cuts of 1.0, 1.5, and 2.0 GeV. The cut on the electron energy is strongly 
correlated with the invariant mass selection on the e/j, pair. Since lower masses 
(in particular m < 4.0 GeV) are thought to have large background contami- 
nation, the loss of these pairs with higher electron energy threshold are not as 
worrisome as they might otherwise be. It may be necessary to use the electron 
energy in order to selectively trigger on these events. 
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Figure 31. DD data from 7r~ — Cu collisions at ^/s = 26 GeV. In particular, panel d) shows 
the DD rapidity difference Ay = yo — yjj- 



Measuring charm by this method requires an accurate model of the cc 
distribution in Q 2 , and the corresponding Ay (rapidity gap) and Ap±. In 
Figure 31 PYTHIA is compared to some of the only data on the rapidity gap 
between D and D mesons [44]. The agreement is not bad, but it the compar- 
ison is not a great confidence builder in the ability to model and then test by 
data checking this input. More experimental data are needed and theoretical 
work to model charm production. 

In order to draw a full picture of charm production, multiple measure- 
ments must be done in the single lepton (electrons and muons) and correlated 
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leptons (electron pairs, muon pairs, and electron-muon pairs). It would be 
extremely useful for theorists interested in total charm or beauty production, 
and also high p± energy loss of heavy flavor partons, to make some predic- 
tions for the leptonic signatures that will be measured in the next year. Future 
upgrades to the detectors for tagging displaced vertices from D meson decays 
are probably more than five years away. 

7.2. QUARKONIA (PHENIX) 




2 4 6 8 10 12 

Muon Pair Mass (GeV) 
Figure 32. Simulation study of dimuons in the PHENIX muon spectrometers. 

There are no early measurement results on quarkonia J/ij) states from 
Run I due to the low luminosity and short running period. The PHENIX ex- 
periment will make a measurement of J/ ip and other states in both the muon 
and electron channels in Run II. PHENIX has a large acceptance in xf and 
p± that we be crucial to constraint models of coloring screening absorbtion 
and test theories with re-coalescence at the hadronization phase. Shown in 
Figure 32 is a simulation of the type of measurement that could be made 
in the PHENIX muons arms with 37 weeks at 10% of design luminosity, or 
equivalently in less than four weeks with the luminosity averaging the the 
RHIC design specification. The first measurements are being made now in 
Run II, and high statistics should be available in Run III. 
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7.3. QUARKONIA (STAR) 



The STAR detector system is unique among the RHIC experiments in its 
capabilities to simultaneously measure many experimental observables on a 
event-by-event basis such as energy-density, entropy, baryochemical poten- 
tial, strangeness content, temperature, and flow. The measurement of J/ty 
production as a function of these quantities allows the study of the sup- 
pression mechanism in great detail. This advantage becomes immediately 
apparent when it comes to the study of the onset of the anomalous suppression 
which reflects the point where the system reaches critical conditions and, 
at least partially, undergoes a phase transition. Thus, the correlation of the 
cc break-up with the many single-event variables provides a new promising 
analysis tool. 

The golden J/* decay mode for STAR is J/ty — > e + e~. Prima facie 
this poses a problem since the STAR detector has been designed to focus 
primarily on hadronic observables over a large phase-space and thus lacks two 
essential features of a dedicated lepton experiment: hadron-blind detectors 
and fast event recording rates. The two essential components which help to 
overcome these shortcomings are ( i) the EMC barrel which allows to suppress 
hadrons to a level sufficient to achieve signal-to-background ratios around 1:3 
or better, and ( ii) a fast level-3 trigger which is designed to efficiently trigger 
on electron-positron pairs with a given invariant mass at rates in the order 
of 100 Hz, thus improving STARs bandwidth for recording J/ty decays by 
almost two orders of magnitude. 

The geometric acceptance for the decay channel J/ty — ► e + e _ in STAR 
is shown in Fig. 33. A J/ty is accepted if both electrons carry momenta 
p > 1.5 GeV/c and fall into the EMC acceptance. Both requirements imply 
that the electron tracks cross all layers of inner tracking detectors (SVT+SSD) 
and all TPC padrows. This ensures maximum momentum and dE/dx resolu- 
tion and therefore maximum additional electron identification from detectors 
other than the EMC. As is depicted in the left plot, the full coverage of the 
EMC \rf\ < 1 ensures relatively high efficiencies (~ 8%) at very low p±. 
At these low values the two leptons run essentially back-to-back, a decay 
topology that requires symmetric coverage around y = 0. This region of 
phasespace is of great interest since here the J/ty remains longest in the hot 
dense medium and its breakup probability is maximal. Up to p± =1 GeV/c 
the acceptance then drops significantly since, still at large opening angles, it 
becomes more likely to loose one electron because it either carries too low 
momentum and/or falls outside our acceptance. At larger p± the acceptance 
raises dramatically due to the decreasing opening angle of the pair an the 
higher average momenta of the electrons. However, the J/* cross-section 
drops exponentially towards larger p± resulting in little net benefit in terms 
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Figure 33. Geometrical acceptance for J/\I> — > e + e _ in the STAR experiment. The left plot 
shows the acceptance as a function of p± for various rapidities y, right one as a function of 
rapidity for various p± slices. 



of total yields. In this region the acceptance scales approximately linear with 
the EMC coverage while at lower p± it scales almost quadratically . 

The dominant background source is the 'combinatorial background' due 
to pions misidentified as electrons. Other sources as electrons from 7r°- and 
7/-Dalitz decays, photo-conversions, and decays of light vector-mesons (p, <fi) 
turned out to be neglectable, since ( i) their average p± is too low and ( ii) their 
rate is small compared to the abundance of misidentified hadrons. In addition, 
most electrons from photo-conversion can easily be rejected by requiring that 
all tracks point back to the primary event vertex. 

The most crucial factor in the overall background rejection is the hadron 
rejection power (e/h) of the EMC. However, for p < 2 GeV/ c, a region where 
the hadron rejection capabilities of the EMC are degraded, the combined 
dE/dx information of SVT and TPC helps to augment the S/B ratio consider- 
ably. It is important to note, that the hadron rejection enters quadratically into 
the background. 

The magnitude of the background depends strongly on the momentum 
cut applied to the electron candidates. In order to access the low-p^ J/* 
region STAR has chosen mj/^/2 ~ 1.5 GeV/c as the lowest value to study. 
The higher the cut the better the e/h rejection from the EMC and the lower 
the charged pion yield. 
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Table I. Current estimates on J/* yields, signal-to-background ratio, and 
statistical significance of the signal in the STAR experiment. 



cut 


J/*-yield 


S/B 


a after 10 7 sec 


a after 10 6 sec 


p e > 1.5 GeV/c 


40 k 


1:3 


76 


24 


p e > 2.0 GeV/c 


10 k 


3:1 


77 


24 



Before any measurement can be addressed it is important to estimate 
the achievable yields and the statistical significance of the signal under re- 
alistic assumptions. Table 7.3 summarizes the resulting yields, the signal- 
to-background ratio, and the statistical significance after 10 7 and 10 6 sec. 
Note, that a nominal RHIC run has 10 7 sec. Here one assumes 100% level-3 
trigger efficiency running at the design rate of 100 Hz, and full EMC cover- 
age including pre-shower detectors. The estimate includes the reconstruction 
efficiencies of the various detectors. Known cross-section from elementary 
collisions were used to exptrapolate to Au+Au at ^/s NN = 200 GeV. Two 
different sets of cuts on the electron momenta are shown in order to demon- 
strate its effect on the S/B ratio. Interestingly, the higher S/B ratio for the 
larger p± cut balances the decreased signal strength, thus resulting in the 
same statistical significance for both cases. 

For the actual measurement one has to study J/ \I/ production not only 
in the most central collisions, but has to vary the centrality and possibly the 
colliding systems. Although the J/* yield will decrease for semi-central 
and peripheral collisions the signal-to-background ratio will decrease even 
stronger (almost quadratically) which will ease the extraction of the signal 
significantly. 



8. Summary and Conclusions 

The workshop was a wonderful forum for learning and exchanging new ideas 
about the physics relevant at RHIC. There is lots of exciting physics in a field 
with much potential for discovery. We wish to thank the organizers and all 
the students for their active participation and thought provoking questions. 
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